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Abstract. In this paper, we derive a general asymptotic implied volatility at 
the first-order for any stochastic volatility model using the heat kernel expansion 
on a Riemann manifold endowed with an Abelian connection. This formula is 
particularly useful for the calibration procedure. As an application, we obtain an 
asymptotic smile for a SABR model with a mean-reversion term, called A-SABR, 
corresponding in our geometric framework to the Poincare hyperbolic plane. 
When the A-SABR model degenerates into the S'A.Bit'-model, we show that our 
asymptotic implied volatility is a better approximation than the classical Hagan- 
al expression |19| . Furthermore, in order to show the strength of this geometric 
framework, we give an exact solution of the SABR model with (3 = or 1. In a 
next paper, we will show how our method can be applied in other contexts such 
as the derivation of an asymptotic implied volatility for a Libor market model 
with a stochastic volatility (|23|). 
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1 Introduction 

Since 1973, the Black-Scholes formula [H] has been extensively used by traders in 
financial markets to price options. However, the original Black-Scholes derivation 
is based on several unrealistic assumptions which are not satisfied under real 
market conditions. For example, in the original Black-Scholes framework, assets 
are assumed to follow log-normal processes (i.e. with a constant volatility). This 
hypothesis can be relaxed by introducing more elaborate models called local and 
stochastic volatility models (see ^2] f° r a n i ce review). 
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On the one hand, local volatility models assume that the volatility <T£,(/, t) de- 
pends only on the underlying / and on the time t. The market is still complete 
and, as shown by Dupire HJ, there is a unique diffusion term (Tl(/, t) which 
can be calibrated to the current market of European option prices. On the 
other hand, stochastic volatility models assume that the volatility itself follows 
a stochastic process [2] : in this case, the market becomes incomplete as it is 
not possible to hedge and trade the volatility with a single underlying asset. It 
can be shown that the local volatility function represents some kind of average 
over all possible instantaneous volatilities in a stochastic volatility model |12| . 

For these two types of models (local and stochastic), the resulting Black-Scholes 
partial differential equation becomes complicated and only a few exact solutions 
are known QJ[2]. The most commonly used solutions are the Constant Elasticity 
of Variance model (CEV) which assumes that the local volatility function is given 
by <Tl(/, t) = (To/ with ao and /3 constant, and the Hcston model which assumes 
a mean-reverting square-root process for the variance. 

In all other cases, analytical solutions are not available and singular perturbation 
techniques have been used to obtain asymptotic expressions for the price of 
European-style options (see ^2 f° r a review of singular perturbation techniques) . 
As for fast reverting volatility processes, similar expansions also exist 1 1 1| . 

These singular perturbation techniques can also be used to derive an implied 
volatility (i.e. smile). By definition, this implied volatility is the value of the 
volatility that when put in the Black-Scholes formula, reproduces the market 
price for a European call option. In the authors discover that the lo- 

cal volatility models predict the wrong behavior for the smile: when the price 
of the underlying decreases (increases), local volatility models predict that the 
smile shifts to higher (lower) prices. This problem can be eliminated with the 
stochastic volatility models such as the SABR model (depending on 4 parame- 
ters: Sigma, Alpha, Beta, Rho) The SABR model has recently been the 
focus of much attention as it provides a simple asymptotic smile for European 
call options, assuming a small volatility. 

Outline 

In this paper, we obtain asymptotic solutions for the conditional probability and 
the implied volatility up to the first-order with any kind of stochastic volatil- 
ity models using the heat kernel expansion on a Ricmannian manifold. This 
asymptotic smile is very useful for calibration purposes. The smile at zero-order 
(no dependence on the expiry date) is connected to the geodesic distance on a 
Riemann surface. This relation between the smile at zero-order and the geodesic 
distance has already been obtained in and used in [2] to compute an asymp- 
totic smile for an equity basket. Starting from this nice geometric result, we show 
how the first-order correction (linear in the expiry date) depends on an Abclian 
connection which is a non-trivial function of the drift processes. 
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We derive the unified asymptotic implied volatility in two steps: 

First, we compute the local volatility function associated to our general stochas- 
tic volatility model. It corresponds to the mean value of the stochastic volatility. 
This expression depends on the conditional probability which satisfies by defi- 
nition a backward Kolmogorov equation. Rewriting this equation in a covariant 
way (i.e. independent of a system of coordinates), we find a general asymptotic 
solution in the short-time limit using the heat kernel expansion on a Rieman- 
nian manifold. Then, an asymptotic local volatility at first-order is obtained 
using a saddle-point method. In this geometric framework, the stochastic (lo- 
cal) volatility model will correspond to the geometry of (real) complex curves 
(i.e. Riemann surfaces). In particular, the SABR model can be associated to 
the geometry of the (hyperbolic) Poincare plane H 2 . This connection between 
H 2 and the SABR model has been obtained in an unpublished paper [27j and 
presented in |28| . Similar results can be found in [S]- 

The second step consists in using a one-to-one asymptotic correspondence be- 
tween a local volatility function and the implied volatility. This relation is de- 
rived using the heat kernel on a time-dependent real line. 

Next, we focus on a specific example and derive an asymptotic implied volatility 
at the first-order for a SABR model with a mean-reversion term which we call 
A-SABR. The computation for the smile at the zero-order is already presented in 
and a similar computation for the implied volatility at the first-order (without 
a mean- reversion term) is done in |2()j . 

Furthermore, in order to show the strength of this geometric approach, we obtain 
two exact solutions for the conditional probability in the SABR model (with 
j3 = ,1). The [3 = solution has already been obtained in an unpublished 
paper |27| and rederived by the author in |22j . For (3 = 1, an extra dimension 
appears and in this case, the SABR model is connected to the three-dimensional 
hyperbolic space H 3 . This extra dimension appears as a (hidden) Kaluza-Klcin 
dimension. 

As a final comment for the reader not familiar with differential geometry, we 
have included a short appendix explaining some key notions such as manifold, 
metric, line bundle and Abelian connection (see |1U| for a quick introduction). 
The saddle-point method is also described in the appendix. 

2 Asymptotic Volatility for Stochastic Volatility 
Models 

A stochastic volatility model depends on two SDEs, one for the asset / and one 
for the volatility a. Let denote x = (xi)i = i.2 = (/, ft), with initial conditions 
a = (ai)j=i,2 = (fo,ct). These variables xt satisfy the following stochastic 
differential equations (SDE) 
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dx l = b t {x,t)dt + <r l (x,t)dW l (2.1) 
dWidWj = pij (x, t)dt (2.2) 

with the initial condition Xi(t = 0) = a, (with b' = in the risk- neutral measure 
as / is a traded asset). 

We have that the local volatility is the mean value of the stochastic volatility 

2, f \ Ia° vi(T,x) 2 p(x,T\a)dx 2 

& {f, T ) = Focr-; rri K 2 -*) 

J p(x,r|a)c(X2 

with p(x,r\a) the conditional probability. 

In order to obtain an asymptotic expression for the local volatility function, 
we will use an asymptotic expansion for the conditional probability p(x,r\a) 
satisfying the Backward Kolmogorov equation 



dp 



£=b i d i p + g i id ijP , (i,j)=a,f (2.4) 



P (t = 0) = S(x - a) (2.5) 

with g % i = I} did j. In (|2.4|l . we have used the Einstein summation conven- 
tion, meaning that two identical indices are summed. For example, dib l p = 
J2i=i 9ib l p. Note that in the relation g 1 ^ = ^Pij^i^j although two indices are 
repeated, there are no implicit summation over i and j. As a result, g lJ is a 
symmetric tensor precisely dependent on these two indices. We will adopt this 
Einstein convention throughout this paper. 

In the next section, we will show how to derive an asymptotic conditional proba- 
bility for any multi-dimensional stochastic volatility models (j2.1J) using the heat 
kernel expansion on a Ricmannian manifold (we will assume here that n is not 
necessary equal to 2 as it is the case for a stochastic volatility model.). In 
particular, we will explain the DeWitt's theorem which gives the asymptotic 
solution to the heat kernel. An extension to the time-dependent heat kernel will 
be also given as this solution is particulary important in finance to include term 
structure. 



3 Heat Kernel Expansion 

3.1 Heat kernel on a Riemannian manifold 

In this section, the partial differential equation (PDE) 12. 4|) will be interpreted 
as the heat kernel on a general smooth n-dimcnsional manifold M (here we 
have that i, j = 1 • • • n) without a boundary, endowed with the metric <?y (see 
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Appendix A for definitions and explanations of the underlined names). The 
inverse of the metric g %3 is defined by 

g 13 = -^Pij^j (3.1) 

and the metric (p y inverse of pij, i.e. p %3 pj k = 8 l k ) 

p tj 

9ij = 2^— (3.2) 

The differential operator 

D = tfdi + g ij dij (3.3) 

which appears in (|2.4(l is a second-order elliptic operator of Laplace type. We 
can then show that there is a unique connection V on C, a line bundle over M, 
and a unique smooth function Q on M such that 

D = ;/"V,V, • () (3.4) 
= g-^{d i +A)9^9 ij {dj+A j ) + Q (3.5) 

Using this connection, (|2.4[1 can be written in the covariant way, i.e. 

d 

g^p(x,y,T) = Dp(x,y,T) (3.6) 

If we take Ai = , Q = then D becomes the Laplacc-Bcltrami operator (or 
Laplacian) A = g~ z dig* g n dj . For this configuration, (|3.6(l will be called the 
Laplacian heat kernel. 

We may express the connection A % and Q as a function of the drift bi and the 
metric gij by identifying in (|3.5|) the terms di and dij with those in (|2.4|l . We 
find 



A* = -(V-g-ifyig 1 / 2 ^)) (3.7) 
Q = ir'iA.A, b,A, iKA,) (3.8) 

Note that the Latin indices i,j ■ ■ ■ can be lowered or raised using the metric gij 
or its inverse g v . For example Ai = gij A 3 . 

The components Ai = gijA 3 define a local one-form A = Aidx 1 . We deduce that 
under a change of coordinates x % (x J ), Ai undergoes the vector transformation 
Ai'diX 1 = Ai- Note that the components bi don't transform as a vector. This 
results from the fact that the SDE (|2.1|) has been derived using the Ito calculus 
and not the Stratanovich one. 

Next, let's introduce the Christoffel's symbol r£j- which depends on the metric 
and its first derivatives 
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2# P ( d j9w + 9i9jp - d p gji) 



(3.9) 



1(3.7(1 can be re-written 



(3 
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Note that if we define 



V = 



e x(^iT)-x(a=oi,r=0) 



P 



(3 
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then p' satisfies the same equation as p ((3.6(1 but with 



A' t = A-da 

Q' = Q+d T X 



(3 



12) 
13) 



The transformation (13.11(1 is called a gauge transformation. The reader should 
be aware that the transformation 1(3.12(1 only applies to the connection Ai with 
lower indices. 

The constant phase e x ^ x=ot ' T=0 ' has been added in 1(3.11(1 so that p and p' satisfy 
the same boundary condition at t = 0. Mathematically, ((3.11(1 means that 
p is a section of the line bundle C and when we apply a (local) Abelian gauge 
transformation, this induces an action on the connection A ((3.12(1 (see Appendix 
A). In particular, if the one-form A is exact, meaning that there exists a smooth 
function \ such that A = d\ then the new connection A' ((3.12(1 vanishes. 

The asymptotic resolution of the heat kernel ((3.6(1 in the short time is an im- 
portant problem in theoretical physics and in mathematics. In physics, it corre- 
sponds to the solution of the Euclidean Schrodinger equation on a fixed space- 
time background [7] and in mathematics, the heat kernel - corresponding to the 
determination of the spectrum of the Laplacian - can give topological informa- 
tion (e.g. the Atiyah-Singer's index theorem) |13] , The following theorem proved 
by Minakshisundaram-Pleijel-De Witt-Gilkey gives the complete asymptotic so- 
lution for the heat kernel on a Ricmannian manifold. 

Theorem 3.1 Let M be a Riemannian manifold without a boundary. Then for 
each x € M , there is a complete asymptotic expansion 



Kx.i/.t) = i^yAMP(„) e -^^ ffln ( Iiy ) T « , r^O (3.14) 



• Here, c(x, y) is the Synge world function equal to one half of the square of 
geodesic distance \x — y\ g between x and y for the metric g. This distance 
is defined as the minimiser of 



oo 



G 



i 1 2 • f dx dx^ , . 

F - V\g = ™ m c J 9ij-^-~^dt (3.15) 

and t parameterises the curve C . The Euler-Lagrange equation gives the 
following geodesic differential equation which depends on the Christoffel's 
coefficients T l jk (|3.9|) 

w + r *^- <"«> 

• A(x, y) is the so-called Van Vleck-Morette determinant 

A(x, y) = g(x)-Uet(- d2 ^ X dy y) )g(y)-i (3.17) 

with g(x) = detgij{x,x) 

• V{x, y) is the parallel transport of the Abelian connection along the 
geodesic from the point y to x 



V(x,y)=e I °(*.v) Aidxi (3.18) 

• The a,i{x,y) (ao(x,y) = 1) are smooth functions on M and depend on 
geometric invariants such as the scalar curvature R. The first coefficients 
are fairly easy to compute by hand. Recently a n {x, x) has been computed 
up to order n = 8. These formulas become exponentially more complicated 
as n increases. For example, the formula for ae has 46 terms. The first 
coefficients are given below 



a x (x,x) = P=-R + Q (3.19) 
6 

a 2 (x,x) = -±-{R l]kl R l3kl - R l0 R lJ + AR 2 ) (3.20) 

loU 

+ \p 2 + ^n rj n^ + ^ap (3.2i) 

with Rijki the Riemann tensor, Rij the Ricci tensor and R the scalar 
curvature given by 

RjM = di^jk - <9 fe r *, + r^Tp/ - r^r pA , (3.22) 

Rjl = R)a (3-23) 

R = g ij Rij (3.24) 

Hij = [Vi.Vj] (3.25) 
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Remark 3.2 In financial mathematics, whenp is computed at n order, we call 
the solution an asymptotic solution of order 5 as the asymptotic order is defined 
according to the dimensionless parameter (a^t) 71 . For example, the SABR model 
(that we will review) is an asymptotic solution at the second-order and only 
requires the knowledge of the first heat kernel coefficient a\ . 

Let's now explain how to use the heat kernel expansion (j3.14|) with a simple 
example, namely a log-normal process. 

Example 3.3 (Log-normal process) The SDE is df = aofdW and using the 
definition ^3.ty) . we obtain the following one- dimensional metric 



9ff = ( 3 ' 26 ) 

When written with the coordinate s = \/2 ln ^ , the metric is flat g ss = 1 and 
all the heat kernel coefficients depending on the Riemann tensor vanish. The 
geodesic distance between two points s and s' is given by the classical Euclidean 
distance d(s, s r ) = \s — s'\ and the Synge function is a(s, s') = ^(s — s') 2 . 
In the old coordinate [f], a is 

*(f,fo) = \ln{^f (3.27) 
CT o fa 

Furthermore, the connection A \3.7]) and the function Q \S. S\) are given by 



A = -Yf d! (3 ' 28) 

Q = -f (3-29) 

Therefore the parallel transport V(f, fo) is given by V(f, fo) = e 2ln ^-Ta\ 
Plugging these expressions into {3.1$ , we obtain the following first-order asymp- 
totic solution for the log-normal process 

/(27rcr^r)2 8 
At the second-order, the second heat kernel coefficient is given by a-i = \Q 2 = 

<T 4 T 2 

° and the asymptotic solution is 



P2(f,r\f ) 



f(2ira 2 T)-2 



128 



(3.31) 
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Furthermore, we note that A and Q are exact, meaning that A = dA, Q = —d T A, 

with A = —^ln(f) + ■ Modulo a gauge transformation p' = n ^7°'~* r p, p' 
satisfies the (Laplacian) heat kernel on IR 



d 2 s p' = d T p (3.32) 

whose solution is the normal distribution. So the exact solution for p = 
-l-ln( t ) °° T 

p g 2 yj^l 8 jg gi ven fry 



p(/,r|/o) = 



1 



f(2na 2 T)- 



(3.33) 



To be thorough, we conclude this section with a brief overview of the derivation 

of the heat kernel expansion. 

We start with the Schwinger-DeWitt antsaz 



P(t,x\v) = ^Ldet[A{x,y)]^e-^V{x,y)^{t,x\y) (3.34) 

Plugging H3.34fl into the heat kernel equation i|3.6[) . we derive a PDE satisfied 
by the function £l(t, x\y) 



d t n = (--a'Vi +V~ 1 A-^(D + Q)A^V)n (3.35) 

with V, = di + Ai and o~i = Vi<7, a 1 = g tJ '<jj. The regular boundary condition 
is Q(t, x\y) = 1. We solve this equation by writing the function 17 as a formal 
series in t: 



oo 

n{t,x\y) = J2 a n(x\y)t n (3.36) 

n=Q 

Plugging this series (|3.36() into (|3.35() and identifying the coefficients in t n , we 
obtain an infinite system of ordinary coupled differential equations: 



ci = f (3.37) 

(f + -jy V 4 H = V^A-i(D + Q)A^Va k -x (3.38) 

The calculation of the heat kernel coefficients in the general case of arbi- 
trary background offers a complex technical problem. The Schwinger-DeWitt 's 
method is quite simple but it is not effective at higher orders. By means of it 
only the two lowest order terms were calculated. For other advanced methods 
see jnj- 



9 



3.2 Heat kernel on a time-dependent Riemannian mani- 
fold 



In most financial models, the volatility diffusion-drift terms can explicitly depend 
on time. In this case, we obtain a time-dependent metric and connection. It is 
therefore useful to generalise the heat kernel expansion from the previous section 
to the case of a time-dependent metric defined by 



This is the purpose of this section. 
The differential operator 

V = b i (t)d i +g ij (t)d ij (3.40) 

which appears in l|2.4J) is a time dependent family of operators of Laplace type. 
Let ; denote the multiple covariant differentiation according to the Levi- 
Civita connection (t). We can expand D in a Taylor series expansion in t to 
write T> invariantly in the form 



Vu = Du + J2 t r Gi j U;ij + (3.41) 



r>0 



with the operator D depending on a connection At and a smooth function Q 
given by l|3.5J) and the tensor Q l J and T % r up to the second-order by 



Q\ 3 = = ^WK-(O) + PiM°iM°M (3-42) 

<?? = \9 l .m - ^f^W^O) + ^1^(0)^.(0) (3.43) 

+ ^^(0K,(0) (3.44) 

T\ = V t (0) (3.45) 

n = b%(0) (3.46) 
Using this connection, (j2.4|l can be written in the covariant way, i.e. 



—p(x,y,t)=Vp{x,y,t) (3.47) 

with V given by (|3~4*T|) . 

The asymptotic resolution of the heat kernel 13.47fl in the short time in a time- 
dependent background is an important problem in quantum cosmology. When 
the spacetime slowly varies, the time-dependent metric describing the cosmolog- 
ical evolution can be expanded in a Taylor series with respect to t. The index 
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r in this situation is related to the adiabatic order [S] . The following expression 
obtained in [141 ITS] gives the complete asymptotic solution for the heat kernel 
on a Riemannian manifold. 

Theorem 3.4 Let M be a Riemannian manifold without a boundary and a time- 
dependent metric. Then for each x € M , there is a complete asymptotic expan- 
sion 



p(x, y, t) = ^^^A(x~y~)P(x, y)e~^ V a n (x, y, t)t n (3.48) 

The a,i(x,y,t) are smooth functions on M and depend on geometric invariants 
such as the scalar curvature R. The coefficients a n have been computed up to 
the fourth-order (ao(x, y, t) = 1). The first coefficients are given below 

1 3 

ai(x,x,t) = P = -R + Q + -Gi,ii 

a 2 (x,x,t) = -L(R ljkl R^ kl - /,',,//" + AR 2 ) + ]-P 2 + -i-7C,,7C'-' + \ap 

1 45 45 
+ -^r^{-^Gi,uGi,jj + —0] .,0] ., + 6OG2.U + 15Si,«-R — : >'^'; ., />',,. 

• (iO.'A'i , : , + l5Qi,M;jj — : >()C/'| i 

4 Geometry of Complex Curves and Asymptotic 
Volatility 

In our geometric framework, a SVM corresponds to a complex curve, also called 
Riemann surfaces. Using the classification of conformal metric on a Ricmann 
surface, we will show that the SVM falls into three classes. In particular, the 
A-SABR and Heston models correspond to the Poincare hyperbolic plane. This 
connection between the SABR model and H 2 has already been presented in 
|51 121)1 1771 128| . This identification allows us to find an exact solution to the SABR 
model (A = 0) with f3 = 0, 1. The (3 = solution has already been obtained in an 
unpublished paper [23121] and rederived by the author in |22J - Furthermore, we 
will derive a general asymptotic implied volatility for any stochastic volatility 
models. This expression only depends on the geometric objects that we have 
introduced in this section (i.e. metric, connection). 

4.1 Complex curves 

On a Riemann surface we can show that the metric can always be written locally 
in a neighborhood of a point (using the right system of coordinates) 



9ij 



«*>ri5ij,i,j=l,2 (4.1) 
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and it is therefore locally conformally flat. The coordinates called the 

isothermal coordinates. Furthermore, two metrics on a Riemann surface, gtj 
and hij (in local coordinates), are called conformally equivalent if there exists a 
function 4>{xi) such that 

9H = e+hy (4.2) 

The following theorem follows from the above observations: 

Theorem 4.1 (Uniformisation) Every metric on a simply connected Rie- 
mann surface 1 is conformally equivalent to a metric of constant scalar curvature 
R: 

1. R = +1: the Riemann sphere S 2 

2. R — 0: the complex plane C 

3. R — — 1: the upper half plane H 2 = {z£ C\Im(z) > 0} 

By the uniformisation theorem, all surfaces falls into these three types and we 
conclude that there are a priori three types of stochastic volatility models (mod- 
ulo the conformal equivalence). In the following, we compute the metric associ- 
ated with the A-SABR model and find the corresponding metric on H 2 In 
this way, the A-SABR model can be viewed as an universal stochastic volatility 
model. Furthermore, we show that the Hcston and A-SABR models belong to 
the same class of conformal equivalence. 

In the next section, we present our general asymptotic implied volatility at the 
first-order and postpone the derivation to subsection 4.3. 

4.2 Unified Asymptotic Implied Volatility 

The general asymptotic implied volatility at the first-order (for any (time- 
independent) stochastic volatility models) , depending implicitly on the metric 
9ij l|3.2f> and the connection Ai (|3.7|) on our Riemann surface, is given by 



„ (Km A) H ^ ft I 9 " T ( 3 M f (c) ,2 , fy"® M jU 

J f<> y/2gff(c) 

(gff'(c)((ln(AgV 2 )'(c) - ^M) + g ff "(c))))) (4.4) 



G 



gff(cW{cy y wu V y ; w 0"(c) 

with c the volatility a which minimizes the geodesic distance d(x,x') on the 
Riemann surface (<fi(x, x') = d 2 (x, x 1 )). A is the VanVleck-Morette determinant 
l|3.17|l . g is the determinant of the metric and V is the parallel gauge transport 
H3.18fl . Here the prime symbol ' indicates a derivative according to a. 

lr The non-simply connected Riemann surfaces can also be classified by taking the double 
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This formula l|4.4|l is particularly useful as we can use it to calibrate rapidly 
any SVM. In the section 5, we will apply it to the A-SABR model. In order to 
use the formula, the only computation needed is the calculation of the geodesic 
distance. For example, for a n-dimcnsional hyperbolic space H n , the geodesic 
distance is known. We will see that H 2 corresponds to a SABR model with 
drift. The general case H n will be used in a future paper [53] to compute an 
asymptotic volatility formula for a SABR-BGM model. 



4.3 Derivation 

4.3.1 Asymptotic probability 

We now have the necessary data to apply the heat kernel expansion and deduce 
the asymptotic formula for the probability density at the first-order for any 
(time-independent) stochastic volatility model. We obtain 

P= 7 T^ T v / A(x, 2 /)P( a; ,y)e- i ^(l + (Q + |)T) (4.5) 

(47TT) o 

We will now derive an asymptotic expression for the implied volatility. The 
computation involves two steps. The first step as illustrated in this section 
consists in computing the local volatility <r(/, r) associated to the stochastic 
SABR model. In the second step (see next section), we will deduce the implied 
volatility from the local volatility using the heat kernel on a time-dependent real 
line. 

We know that the local volatility associated to a SVM is given by 



Jo V^fgda 

with p the conditional probability given in the short time-limit at the first-order 
by (IQll . 

We set <j)(z,z') = d 2 (z,z'). The factor ^fg is the (invariant) measure (see Ap- 
pendix A) and g^ = with the definition (13.211 . 

a f 

Plugging our asymptotic expression for the conditional probability (|4.5|l in Ij4.6|l . 

we obtain 



~ / oo ft(o,r)e»*(«)da 1 ' 



with 4>{a) = d 2 (z,z'), h{a,r) = ^g^/A(z, z')V(z, z'){\ + (Q + f )r), /(a,r) = 
h(a,r)g ff and x = -i. 

Using a saddle-point method, we can find an asymptotic expression for the lo- 
cal volatility. For example, at the zero-order, a 2 is given by g^(c) with c the 
stochastic volatility which minimises the geodesic distance on our Riemann sur- 
face: 
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c = a | min a (j) (4-8) 

Using the saddle-point method at the first-order, we find the following expression 
for the numerator in (|4.7|) (sec Appendix for a sketch of the proof) 



r h \ ^Wrf = / 2n f ( vi + ! r /"( c ) ^ (4) (c) fjewie) _ nr(c)) 2 „ 

Jo [ ' y-^"(c) An z l 2/(c)0"(c) + 80"(c)2 + 20»(c)2/(c) 24(</>"(c))3 Jj 

Computing the denominator in i|4.7[l in a similar way, we obtain a first-order 
correction of the local volatility 



U ' j ~ 5 (CK1+ x ( 2<P"(c) { f(c) h(c) )+ 2r(c)i [ f(c) h{c) ))} 
Plugging the expression for / and g, we finally obtain 



- V^(l + ^(^| ((i „(A^)'(„, - + ^|)))(4.9) 

Here the prime symbol ' indicates a derivative according to a. This expression 
depends only on the metric and the connection A on our Riemann surface. 
The final step is to obtain a relation between a local volatility function tr(/, t) 
and the implied smile. We will show in the next section how to obtain such a 
relation using the heat kernel expansion on a time-dependent one-dimensional 
real line 

4.3.2 Local Volatility Model and Implied Volatility 

Let's assume we have a local volatility model 



df = C(f,t)dW t f(t = 0) = f (4.10) 

The fair value of a European call option (with maturity date r and strike K) is 
given by (using an integration by parts and assuming that r is small) 



C(K, r, / ) = (/o - K)+ + ° 2{K ^ 0) ^ dTp(K, T\f Q ) (4.11) 
with p(K,T\fo) the conditional probability. 

In our framework, this model corresponds to a (one-dimensional) real curve 
endowed with the time-dependent metric gff = c(j t yi ■ F° r t = and for 
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the new coordinate u = V2 f ^jt\ (with C(/, 0) = C(f)), the metric is flat: 
g-uu = 1- The distance is then given by the classical Euclidean distance 

d(u,u') = \u- u'\ (4.12) 

Furthermore, the connection A <|3.7|) and the function Q 1)3. 8f) are given by 

A, = -\d f ln{C{f)) (4.13) 

Q = ^&- 1 2& ] (4 ' 14) 



The parallel transport is then given by V(f, fo) = y c{j) 
Furthermore, Q is given by l)3.42[l 



G(K, T) = 2d t ln(C(K, t))\ t=0 (4.15) 

The first-order conditional probability (using the heat kernel expansion on a 
time-dependent manifold l|3.48|l is then 



Plugging this expression in 1)4.11(1 . the integration over t can be performed and 
we obtain 

Proposition 4.2 

C(K,tJo) = (fo -K)+ + ^ ( fg /o)T ( gl M + (Q(K) + 3 -^l)rH 2 (^.m) 



fTi(u) = 2(e-^ 2 + Vm^(N(y/2\u\) - 1)) (4.17) 

(1- 

df 



H 2 (uj) = l(e- uj2 (l~2Lj 2 )-2\uj\ 3 ^(N(V2\uj\)-l)) (4.18) 



! fo V*rC(f>) 

In the case of a constant volatility, the above formula reduces to 
Example 4.3 (Black-Scholes Vanilla option) 



(4.19) 



C(K,r,fo) = (f Q K)+ + y^Ml {HliQ) + {Q{K) + M^) )Tjff2( ^ ))( 4.20) 

2v 2ir 4 

with Oj = 

V 2tcto 
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By identifying the formula (|4.16|> with the same formula obtained with an implied 
volatility a = a{K,T) (|4.20() . we deduce 



At the zero-order, we obtain lu = uj i.e. 



o% s (.K,T) = - g ^ r - (4.22) 

Jfo C(f') 

The formula l|4.22[) has already been found in 0] and we will call it the BBF 
relation in the following. Then using the recurrence equation Ij4.21|l . we obtain 
at the first-order 



°BS{K,T) = -^-(1 + \{\{-JT^f + W + ^f 1 ) (4-23) 

Jfo C(f') Jfo C(f') 

l <fo ) n , C\f)T C"{f) ( C'(f) 2 1 d t C 

Jfo C(f') 

In the case C(f,t) = C(f), we reproduce the asymptotic implied volatility ob- 
tained by Hagan- Woodward (EH)- Now plugging the local volatility l|4.9|) into 
the implied volatility l|4.24|) , we find l|4.4|) and this achieves our derivation of an 
asymptotic implied volatility at the first-order. 



5 A-SABR model and hyperbolic geometry 
5.1 A-SABR model 

The volatility a is not a tradable asset. Therefore, in the risk-neutral measure, 
a can have a drift. A popular choice is to make the volatility process mean- 
reverting. Therefore, we introduce the A-SABR model defined by the following 
SDE 



df = aC{f)dW l 
da = X(a — X)dt + vadWz 
C(f) = f, o(0) = a,/(0) = / 

where W\ and Wi are two Brownian processes with correlation p £ ] — 1, 1[. The 
stochastic Black volatility is a t = af 13 ^ 1 . In the following section, we present 
our asymptotic smile for the A-SABR model and postpone the derivation to the 
next section. 



16 



5.2 Asymptotic smile for the A-SABR 

The asymptotic smile (with strike / , maturity date r and spot f ) at the first- 
order associated to the stochastic A-SABR model is 



^(/o,/,r) = ^(l + «l(^)r) (5.1) 



with 



,,s (aminC(f)) 2 1 2d ff (C{f)a mm ) d f (C(f)a min ) 2 

24 [ p + C(f)a mm 1 C(f)a mm ' ' 

av 2 Tln(Py(a min )(l - p 2 ) y/ cosh(d(a min )) 2 - 1 

with g = (i-ff) (g 7^ 1), « o*(g) - £log( ^i^j ) and 

a TO in(9) = v /q;2 + 2q;i/ P < 3 , + J /2 9 2 - 
Moreover we have 

■p \ i / 

^i -pSABs) ( a ™«0 = ~(Go(#2(fijTim), ^o( a min)i B)0 2 (a m i n ) — Gq (^1 (djnin ) j 4(»min)) -5)^1 



1/ 



+^(a mm )(Gi(0 2 (a mm ),S)-G , i(0 1 (a mm ),B))) (5.3) 

with 



!n(P SAM )'(a min ) = — ^- ^-(Fo(0 2 (a mm ),A( amm ),B)^(a mm ) (5.4) 

2(1 - p z ){\ - p) 

), B)8 / 1 (a min ) - A'(a)(F 1 (0 2 ( ),i?)-F 1 (0 1 ( ),A(a mm ),(B)§) 

and with 



IT 



Gi(x, b) = - csc(x) + 6i?e(log(tan(|))) 



Go(x, a, b) = cot(x) csc(x) (a + sin(x)) (1 + 6 tan(x)) 



=(/) J 

a,, s _ j y/1 -p 2 (ttp+^(/)) 



F (a;,a, 6) = 



c(/r (p (a-c(f)) + vq(f)) 
sin(x) 



a + cos(:r) + b sin(x) 

a, 6) = 



2 6 arctan( fc+ ^±g^ 2) ) -1 + q 2 + a 6 sin(a:) 



(_1 _|_ a 2 _ £2)1 ( — 1 + a 2 — 6 2 ) (a + cos(x) + b sin(x)) 



fT~~Z_ 2 

h(a m in) = 7r - arctan( ) 



?i(a mm ) = -arctan( ■ — ) + rh ap+vq ff))> 

q/ / \ _ <*9i(o m in) 



a ( v l(f) + P (a + a mm )) 




is^T^7q(f) (2ap+vq(f)) a. 
(fa -fa 13 q(f) 



B 



f o (-1 + /3) a mm 

/o!/ ( ap + ^ O (/)) + / /3 a (-1 + ^) (a + ^pg(/)) 
Jo' 3 (-1 + 13) a? min + p (a - a mm )) 

P 



5.3 Derivation 

In order to use our general formula for the implied volatility, we will compute 
the metric and the connection associated to the A-SABR model in the next 
subsection. We will show that the A-SABR metric is diffeomorphic equivalent 
to the metric on H 2 , the hyperbolic Poincare plane [271 128j . 

5.3.1 Hyperbolic Poincare plane 

The metric associated to the A-SABR model is (using (|3.2(l ) 
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ds 2 = gijdx l dx^ (5-6) 

= 2n2 L ^W 2 df 2 -2vpC(f)dadf + C(f) 2 da 2 } (5.7) 

a i (J z v L yY — p z ) 

Let's introduce the variable q(f) = jJ o (Wjn ■ By introducing the new coordi- 
nates x = vq — pa and y = (1 — p 2 )ia, the metric becomes (after some algebraic 
manipulations) the standard hyperbolic metric on the Poincarc half-plane H 2 in 
the coordinates [a;, y] 2 



2 2 dx 2 + dy 2 

ds = — — 2 — ( 5 - 8 ) 

The unusual factor ^ in front of the metric l|5.8[) can be eliminated by scaling 

the time r' = -^-r in the heat kernel (|3.6fl (and Q becomes -jiQ)- This is what 
we will use in the following. 

Remark 5.1 (Heston model) The Heston model is a stochastic volatility 
model given by the following SDEs \21f : 



df = afdWx (5.9) 

da = -Cl + ^ { l- { «f)dt + V -dW 2 (5.10) 
8a 2 a 2 

Let's introduce the variable x = ^ln(f) — p\ , y = (1 — p 2 )^ \ ■ Then, in the 
coordinates [x, y], the metric becomes 



ds2 = -77, ^TTVdsw (5.11) 

?7 (1 — p z ) 2 

and therefore belongs to the same class of conformal equivalence as H 2 . Mod- 
ulo this equivalence, the X-SABR and Heston models correspond to the same 
geometry. 

As this connection between the A-SABR model and H 2 is quite intriguing, we 
investigate some of the useful properties of the hyperbolic space (for example 
the geodesies). First, by introducing the complex variable z = x + iy, the metric 
becomes 



2 How can we prove that this is the correct metric on H 2 ? By applying a Moebius transforma- 
tion (see above), the upper half plane is mapped to the Poincare disk T> = {z S C \z\ < 1}. 

Then if we define x\ = 1 x i = yzy~p"> X3 = 1— fzp ' we °bt am that T> is mapped 

to the Minskowski pseudo-sphere — + x\ + £?, = 1. On this space, we have the metric 
ds 2 = — dxQ + dx\ + cix?,. We can then deduce the induced metric on the Minskowki model. 
On the upper-half plane, this gives 15.81 (without the scale factor -%). 
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IV1 oeta i i_i s Transformation 




Upper half plane 



Figure 1: Poincare disk T> and upper half-plane H with some geodesies. In the 
upper half-plane, the geodesies correspond to vertical lines and to semi circles 
centered on the horizon Im(z) = and in T> the geodesies are circles orthogonal 
to V 



dzdz 
Im(z)' 2 



In this coordinate system, it can be shown that PSX(2, R) 
meaning that the distance is preserved. The action of PSL(2, 
tive and given by 



(5.12) 

is an isomctry, 
'.) on z is transi- 



cz + d 



(5.13) 



Furthermore, let's define the Moebius transformation T as an element of 
PSL(2,M.) which is uniquely given by its values at 3 points: T(0) = 1, T(i) = 
and T(oo) = — 1. If Im(z) > then |T(z)| < 1 so T maps the upper half-plane 
on the Poincare disk V = {z £ C | \z\ < 1}. In the upper half-plane, the 
geodesies correspond to vertical lines and to semi circles centered on the horizon 
Im(z) = 0, and in T> the geodesies are circles orthogonal to T> (Fig. 



3 PSL(2,R) = 5L(2,R)/Z 2 with 5L(2,R) the group of two by two real matrices with 
determinant one. 1^ acts on A G 5L(2,R) by = —A 
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Solving the Euler-Lagrange equation l|3.16|) . it can be proven that the (hyper- 
bolic) distance (invariant under PSL(2,M)) between the points z = x + iy, 
z' = x' + iy' on H 2 is given by (see Appendix A) 



d(z,z') =cosh~ 1 [l+ |Z Z '} 2 ] (5.14) 

*yy 

Using this explicit expression for the hyperbolic distance, the VanVleck-Morette 
determinant is 



A(z z') = d (z,z') /g 15 % 

y/cOsh 2 {d(z,z')) - 1 

Using our specific expression for the geodesic distance (15.1411 , wc find that a m i n , 
the saddle-point, is the following expression 



a mm {f) = y/a 2 + 2ais pq + i; 2 q 2 (5.16) 

As A depends only on the geodesic distance d and d is minimised for a = a m i n7 

we have Zn(A)'(a m i n ) = 0. 

Furthermore, we have (using Mathematica) 



iii i \ 2d(a m i n ) 17 \ 

a-minty {U>min) — —~ ~ / TTTi N\9 "T (0.1 I ) 

a(l - p 2 )y/cosh{d{a m i n )Y - 1 



4>"'{a min ) 3 



(5.18) 



fl(amm) = COS/l ( — j- ) ( 5 ' 19 ) 



_ lf -qvp - ap 2 + a m m {q) ■ 
a(l-p 2 ) 



This above formula (|5.19(1 has already been obtained in |S] . 

5.3.2 Connection for the A-SABR 

In the coordinates [a, /], the connection A is 



1 , (2\\p-2\pa-va 2 C'(f)) (-2AA + 2Aa + vpa 2 C\f)) , . 

^ = ^ 2T( TvTw H ~2 da > 

2(1 — fr) vC{j)a z 



In the new coordinates [x, y], the above connection is given by 



^ = ^Bfl + ^_v_ ^» 1{pdx _ ^i- p 2 dy) (5.20) 
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with 



P dx 



A = -oTi 2ui m ^ for /3 ^ 1 (5.21) 

2(1 - p 2 )(l -/?)/ , e w + i^B_Jl) 

^ MM = ~w^» s " 1 (5 - 22) 

The pullback of the connection on a geodesic C satisfying [x — Xo(a, f)) 2 + y 2 = 
R 2 (aJ) is given by {{5 £ 1 4 ) 



iM = f/^ + MzMEg±^) ( W + 1)4,(5.23) 

and with 

■»^SASi? = - ydy 24) 

2(1 - P 2)(l - /?) ^R2 (ai /} _ /} + ^ R 2 {a] f)- y 2 + ^^yj' 

with i : C — > H 2 the embedding of the geodesic C on the Poincare plane and 
%o = x a + Twn- We have used that i*dx 



Note that the two constants xo and i? are determined by using the fact that the 
two points z\ = —pa + i\J\ — p 2 a and z% = v jj Q (jfji\ — pa + i\J\ — p 2 a pass 
through the geodesic curves. The algebraic equations given R and xq can be 
exactly solved: 



MaJ)= Xf 9 f 2+y \ Vl (5-25) 
2(xi - x 2 ) 



ma, /) = * fe - ^ 2 + ^ 2 + - - + y' 2 ^ 

2 y (xi - x 2 y 

Using polar coordinates x — xq = Rcos(9) 1 y = Rsin(9), we obtain that the 
parallel gauge transport is 



ln{V){a) = ln(V SABR )(a) + ^- f ^ (1 + Btan(0))( S tn(x) + A ) 

f 7ei(o) szn(x) 2 

with 

-pSABRt P sin(9)d6 



M = ^1 2(1-^(1-/3) (cosW + f + ^.s^ 8) 



The case /3 = 1 will be treated in the next section 
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with 0,(0,/) = arctani^) , i = 1,2, A(aJ) = B = -jf= 2 and 

= ~ P ~ • ^he ^ W0 integration bounds #1 and 82 explicitly depend on a 
and the coefficient ^a. Doing the integration, we obtain i|5.3[) . 
Plugging all these results 1|5. 1715. 1815. 1915. 5)l in l|4.4|l . we obtain our final ex- 
pression for the asymptotic smile at the first-order associated to the stochastic 
A-SABR model (|5~Tj) . 

Remark 5.2 (SABR original formula) We can now see how the classical 
Hagan-al asymptotic smile formula can be obtained in the case A = and 
show that our formula gives a better approximation. 
First we approximate a m i n by the following expression 

ttmm - « + qpv (5.29) 

In the same way, we have 



V 1 cosh(d{a mm )) 2 -_L ^ 1 ^ 



Furthermore, for A = ; the connection L5.20}) reduces to 



A = y—^ {-dln{C{f)) + P -d s Cda) (5.31) 

Therefore, the parallel gauge transport is obtained by integrating this one-form 
along a geodesic C 



V = exp( 1 (-ln(^) + f ?d f Cda)) (5.32) 
2(1 - p 2 ) C(/ ) Jc v 

The component f of the connection is an exact form and therefore has easily been 
integrated. The result doesn't depend on the geodesic but only on the endpoints. 
However, this is not the case for the component A a - But by approximating 
d'fC(f') ~ dfC(f) , the component A a becomes an exact form and can therefore 
be integrated 

J ^-d f Cda ~ ?-d f C(f){a - a) (5.33) 

Finally, plugging these approximations into our formula \4-4\j t we reproduce the 
Hagan-al original formula llftf 

^(/o,/,r)^(l + ^(^)r) (5.34) 
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with 

- —24- { P + C{f) ~ { ~CUT ] ] 4 + _ 2^ 1 ' 

Therefore, the Hagan-al classical formula corresponds to the approximation of 
the Abelian connection by an exact form. The latter can be integrated outside 
the parametrisation of the geodesic curves. 

Remark 5.3 (H 2 -model) In the previous section, we have seen that the X- 
SABR mode corresponds to the geometry ofW 2 . This space is particularly nice 
in the sense that the geodesic distance and the geodesic curves are known. A 
similar result holds if we assume that C(f) is a general function (C(f) = f@ for 
X-SABR). In the following, we will try to fix this arbitrary function in order to 
fit the short-term smile. In this case, we can use our unified asymptotic smile 
formula at the zero-order. The short-term smile will be automatically calibrated 
by construction if 



o-ioc(f) = C(f)a mm (f) (5.35) 
a mm (f) = a 2 + 2pavq + v 2 q 2 (5.36) 

q= Ilw) (5 - 37) 

with o~i oc (f) the short-term local volatility. By short-term, we mean a maturity 
date less than 5 years. Solving 1^5.35)) according to q, we obtain 



uq = -pay+Ja^{-l + (P) + ^l^- (5.38) 



and if we derive under f , we have ( if>(f) = ^ff^P ) 

dip v 



y/ijj 2 -a 2 {l- p 2 ) o-ioc(f) 
Solving this ODE, we obtain that C (/) is fixed to 



df (5.39) 



C{f) = ^ (5.40) 

a^l- p 2 cosh\v J fo ^jTj] 



Using the BBF formula, we have 



C(f) = ^l 9 ™ (5-41) 

ay/1 - p 2 cosh[v aB Jj> f) ] 

Using this function for the X-SABR model, the short term smile is then auto- 
matically calibrated. 
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T=5 y, lambda=-10% lambdabar=0% beta=30% nu=27.48% 
rho=11.48% alpha=1.23%, f=100% 




92.00% 



97.00% 



102.00% 



1 07.00% 



Lambda=-10%, lambdabar=0%, f0=5.37%, tau=5 y, nu=27.48%, 
rho=1 1 .48%, beta=30% 
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Lambda=-10%, Lambdabar=2%, f0=5.37%, beta=30%, 
rho=1 1 .48%, Nu=27.48%, T=5 y 

35 




12% 



Figure 2: Pdf p(K,T\f ) = §t^- Asymptotic solution (pink) vs numerical 
solution (blue) (pde solver). In some cases, the Hagan original formula has been 
plotted to see the impact of the mean-reversion term. 

6 Analytical solution for the SABR model with 

= 0,1 

6.1 SABR model with = and H 2 

For the SABR model, the connection A and the function Q are given by 

A = on 1 2 A -d f ln(C)df + ^d s Cda) (6.1) 
2(1 — p z ) v 

For (3 = 0, the function Q and the potential A vanish. Then p satisfies a heat 
kernel where the differential operator D reduces to the Laplacian on H 2 : 



— = A B 2 P (6.3) 

= y 2 {dl + d 2 y )p (6.4) 

Therefore solving the Kolmogorov equation for the SABR model with j3 = 
(called SARO model) is equivalent to solving this (Laplacian) heat kernel on H 2 . 
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Figure 3: Conditional probability for the SABR model with /3 = versus nu- 
merical PDE. 



Surprisingly, there is an analytical solution for the heat kernel on H 2 found 
by McKean [35]. It is connected to the Selberg trace formula |17j . 
The exact conditional probability density p depends on the hyperbolic distance 
d(z, z') and is given by (with r' = iL ^ L ) 



p{d,r') = 2 2 2-r' ?e * / — r db 



Id(z.z') (coshb — coshd(z, z'))i 
The conditional probability in the old coordinates [a, /] is 

dfda „ s 3,3 t' f°° be 



/n ,„ , v dfda „ 5 3,3 t' /" 
P (f,a,T')dfda=- - T ^_ 2- = 7r-?r'- = e-T / 



d(z,z') (coshb — coshd(z 

We have compared this exact solution (Fig. |3J with a numerical PDE solution of 
the SARO model and found agreement. A similar result was obtained previously 
in |27l for the conditional probability. 

The value of a European option is (after an integration by parts) 

C(f, K) = (f-K)+ + - dr daa 2 p{ Xl = K, a, t\o) 
1 Jo Jo 
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In order to integrate over a we use a small trick: we interchange the order of 
integration over b and a. The half space b > d with a arbitrary is then mapped 
to the half-strip a m i n < a < a max and b > l m i n where 5 



0>max 0>rain 



ly/ (a(ch(b)+ p 2 ~ch{b)p 2 )-vpqf - {a 2 - 2 a v p q + v 2 q 2 ) 



a-rmn + a max = 2ach(b) + 2ap 2 (l - ch(b)) + 2vpq 
num 

cosh{l m in) = -5 — 
den 

with 



num = —a 2 — 2 a 2 p 2 — 2 a 2 p 4 + 4 a 2 p e — 6 a v p 3 s + 4 a v p 5 s — 3 v 2 p s 2 

+2 v 2 p i s 2 + p (-1 + 2 p 2 ) {ap + v s) v /a 2 +2aup (-l + 2p 2 ) s + v 2 (-\ + 2p 2 ) 



den = a (-1 + p 2 ) (l + 2p 2 ) (2p(ap + vs)+ ^/a 2 + 2avp (-l + 2p 2 ) s + v 2 (-1 + 2 p 2 ) s 2 ) 

Performing the integration according to a leads to an exact solution for the 
SABR model with (3 = 



dtr'ie-i I b ^ max amin) ( =db 



v(\-p 2 )i Jo Ji min \J (coshb - coshl min ) 

6.2 SABR model with (3 = 1 and the three-dimensional 
hyperbolic space H 3 

A similar computation can be carried out for j3 = 1. Using l|6.1|l . we can show 
that the potential A is exact, meaning there exists a smooth function A such 
that A — dA with A = 2 (i- P ^) (~^ n (f) + f a )- Furthermore, using (|6.2|) . we have 

n — a 2 - 

_ 8(l-p a ) _ 8(l-p 2 ) 2 ■ 

Applying an Abelian gauge transformation p' = e A p (I3.11|) . we find that p' 
satisfies the following equation 



y 2 ( d " + d l Av 2 {l _ p 2y )P' = d T'P' (6-5) 

How do we solve this equation? It turns out that the solution corresponds 
in some fancy way to the solution of the (Laplacian) heat kernel on the three 
dimensional hyperbolic space H 3 . This space can be represented as the upper- 
half space H 3 = {x = (xi, X2, Xs)\x3 > 0}. In these coordinates, the metric takes 
the following form 

5 all these algebraic computations have been done with Mathematica 
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and the geodesic distance between two points x and x' in H 3 is given by 6 

cosh{d(x,x')) = 1+ \ x ~ x '} 2 ( 6 . 7 ) 
2x 3 x 3 

As in H 2 , the geodesies are straight vertical lines or semi-circles orthogonal to 
the boundary of the upper-half space. An interesting property, useful to solve 
the heat kernel, is that the group of isometrics of H 3 is PSL(2, C) 7 . If we 
represent a point p £ H 3 as a quaternion 8 whose fourth components equal zero, 
then the action of an element g £ PSL(2, C) on H 3 can be described by the 
formula 

' ap + b 

V = 9-P = — 6.8) 

cp + d 

with p = xil + x 2 i + x 3 j. 

The Laplacian on H 3 in the coordinates [xi, Xi, x%\ is given by 



A H3 =xl(d 2 Xi +di 2 +d 2 X3 ) (6.9) 
and the (Laplacian) heat kernel is 

d r ,p' = A n3 p' (6.10) 

The exact solution for the conditional probability density p'(d(x 1 x') 1 1), depend- 
ing on the geodesic distance d(x,x'), is |16| 

(47rr')2 sinh(d(x,x')) 

Let's apply a Fourier transformation on p along the coordinate x\ (or equiva- 
lent ly X2) 

p'(x 1 ,x 2 ,x 3 ,x',t') = - 7 =e lkxi p(k,X2,x 3 ,x',t) (6.12) 

J-00 V27T 



Then p' satisfies the following PDE 

6 1 ■ I is the Euclidean distance in R 3 

7 PS X(2,C) is identical to PSL(2,M.), except that the real field is replaced by the complex 
field. 

8 The quaternionic field is generated by the unit element 1 and the basis i, j, k which satisfy 
the multiplication table i.j = k and the other cyclic products. 
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d T ,p> = xU-k 2 + d 2 X2 +dl)p' (6.13) 

By comparing (j6.13|) with l|fi.5[l . we deduce that the exact solution for the con- 
ditional probability for the SABR model with j3 = 1 is (with x = x%, y = £3, 

^ = 2iy(l-p2) I ^1 = 0) 



I I is ^-(In(^)--S-(a-a)) f°° dX\ 1 d(x,x') 

J-oo V27T (47rr')2 sh{d(x,x')j 

A previous solution for the SABR model with /? = 1 was obtained by 
although only in terms of Gauss hypergeometric series. 



7 Conclusions and future work 

Let's summarize our findings. By using the heat kernel expansion, we have 
explained how to obtain a general asymptotic smile formula at the first-order. 
As an application, we have derived the smile formula for a SABR model with 
a mean-reversion term. Furthermore, we have shown how to reproduce [181 
I19| . In the case of the SABR model with (3 = 0,1, exact solutions have been 
found, corresponding to the geometries of H 2 and H 3 , respectively. We note that 
these solutions are not easy to obtain without exploiting this connection with 
hyperbolic geometry. 

This geometric framework allows to obtain other analytical solutions for stochas- 
tic volatility models. For example, a solution for the heat kernel on a multi- 
dimensional hyperbolic space HP exists. In a future work, we will show that this 
space corresponds to a BGM model coupled with a SABR model. Using the 
general heat kernel expansion on a manifold, we will derive an asymptotic smile 
for this particular model 
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8 Appendix: Notions in differential geometry 
Riemanian manifold 



MANIFOLD 




A real n-dimcnsional manifold is a space which looks like K™ around each 
point. More precisely, M is covered by open sets Ui (topological space) which are 
homcomorphic to l n meaning that there is a continuous map 4>i (and its inverse) 
from Ui to W 1 for each i. Furthermore, we impose that the map fa j = 4>~ 1 0(j)j 
from R" to M" is C°°(R n ). 

As an example, a two-sphere S 2 can be covered with two patches: Um and 
Us, defined respectively as S 2 minus the north pole, and the south pole. We 
obtain the map <f>^ (4>s) by doing a stercographic projection on Un (Ms)- This 
projection consists in taking the intersection of a line passing through the North 
(South) pole and a point p on S 2 with the equatorial plane. We can show that 
4>sn(x, y) = ( u'i^y'is , — (giij^ii) ) is C°° and even holomorphic. So S 2 is a complex 
manifold. 



N 




S 

SPHERE 



Metric 

A metric gij written with the local coordinates Xi (corresponding to a particular 
chart Ui) allows us to measure the distance between infinitesimally nearby points 
x % and x l + dx l : ds 2 = gijdx % dxK If a point p belongs to two charts then the 
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distance can be computed using two different systems of coordinates x % and 
x l = f(x l ). However, the result of the measure should be the same, meaning 
that gijdx l dxi = gijdx 1 dx' . We deduce that under a change of coordinates, the 
metric is not invariant but changes in a contravariant way by c/y = g^jidiX 1 djX J . 
A manifold endowed with an Euclidean metric is called a Riemannian manifold. 
On a n-dimensional Riemannian manifold, the measure y/g Yi7=i dx l is invariant 
under an arbitrary change of coordinates. Indeed the metric changes as gtj = 

gi'j'^r^r and therefore g = det(gij) changes as ^fg ~ det{^-r)\/g 1 . Further- 
more, the element Y\a=i changes as ll"=i = ^ H^t) ll™=i ^ x% an d 
we deduce the result. 

Line bundle and Abelian connection 
Line bundle 



FS 




Line Bundle 



A line bundle C is defined by an open covering of M, {U a }, and for each (a, (3), 
a smooth transition function g a p : U a fl Up — > R which satisfies the " cocycle 
condition" 

9apgf3 1 9 ia = 1 on UaDUpDU-f (8.1) 
A section a of C is defined by its local representatives a on each U a : 

a\ Ua = o~ a '. U a > M (8.2) 
and they arc related to each other by the formula a a = g a p<Jp on U a D Up. 

Abelian connection 

An Abelian connection V on the line bundle £ is a collection of differential 
operator di + Af on each open set U a which transforms according to 

Af = A? + g a pd ig ~l on U a n Up (8.3) 
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9 Appendix: Laplace integrals in many dimen- 
sions 



Let u be a bounded domain on R™, 5 : f2 — > R, / : f2 — > R and A > be a large 
positive parameter. The Laplace's method consists in studying the asymptotic 
as A — > oo of the multi-dimensional Laplace integrals 

F(X) = / f{x)e xs{x) dx (9.1) 
Jn 

Let 5* and / be smooth functions and the function S have a maximum only at 
one interior nondcgcncrate critial point xq £ Q. Then in the neighborhood at 
xq the function S has the following Taylor expansion 

S(x) = S(x ) + \{x- x ) 2 d 2 x S(x )(x - xo) + o((x - x Q ) 3 ) (9.2) 

As A — > oo, the main contribution of the integral comes from the neighborhood 
of xq. Replacing the function / by its value at xq, we obtain a Gaussian integrals 
where the integration over x can be performed. One gets the leading asymptotic 
of the integral as A — > oo 

F(X) ~ e xs ^\^[-det(d 2 x S(x ))}^f(x ) (9.3) 

More generally, doing a Taylor expansion at the n th order for S (resp. n— 2-order 
for f(x)) around x = Xo, we obtain 

2tt -i °° 

F(X) ~ e^(«o)( )i HM(^S(xo))]^5> fe A- fc (9.4) 

fe=0 

with the coefficients ak are expressed in terms of the derivatives of the functions 
/ and S at the point xq. For example, at the first-order (in one dimension), we 
find 



F(X\ I 2n pA5(xo) ff / r ^ , h Ago) i /(^ (4) M , f'(xo)SW(x ) 5S'"(x )f(x ) 

F{X) V ^XS^f {f{Xo) + + 8S»(x )2 + 2S"(x r 245"(x )3 )l(9 ' 5) 
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